Two-dimensional elastic turbulence 
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We investigate the effect of polymer additives on a two-dimensional Kolmogorov flow at very low 
Reynolds numbers by direct numerical simulations of the Oldroyd-B viscoelastic model. We find 
that above the elastic instability threshold the flow develops the elastic turbulence regime recently 
observed in experiments. We observe that both the turbulent drag and the Lyapunov exponent 
increase with Weissenberg, indicating the presence of a disordered, turbulent-like mixing flow. The 
energy spectrum develops a power-law scaling range with an exponent close to the experimental and 
theoretical expectations. 



One of the most remarkable effects of highly viscous 
polymer solutions which has been recently observed in 
experiments is the development of an "elastic turbulence" 
regime in the limit of strong elasticity [![. The flow 
of polymer solution in this regime displays irregularities 
typical of turbulent flows (broad range of active scales 
and growth of flow resistance) even at low velocity and 
high viscosity, i.e. in the limit of vanishing Reynolds 
number. As a consequence of turbulent motion at small 
scales, elastic turbulence has been proposed as an effi- 
cient technique for mixing in very low Reynolds flows, 
such as in microchannel flows 0, 0, 0| ■ Despite its great 
technological interest, elastic turbulence is still only par- 
tially understood from a theoretical point of view. Recent 
theoretical predictions are based on simplified versions of 
viscoelastic models and on the analogy with MHD equa- 
tions 0,0. 

In this letter we investigate the phenomenology of elas- 
tic turbulence in direct numerical simulation of polymer 
solutions in two dimensions. Our main objective is to 
show that usual viscoelastic models, developed for study- 
ing high Reynolds turbulent flows, are able to capture, 
in the limit of vanishing Reynolds numbers, the main 
phenomenology of elastic turbulence, i.e. irregular tem- 
poral behavior and spatially disordered flow. Despite the 
important geometrical differences, our numerical results 
are in remarkable agreement with experimental observa- 
tions of elastic turbulence: this suggests the possibility 
to understand elastic turbulence on the basis of known 
viscoelastic models. 

To describe the dynamics of a dilute polymer solutions 
we adopt the well known linear Oldroyd-B model Q 

ft« + (u ■ V)m = - Vp + vAm + — V ■ u + / (1) 
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where u is the incompressible velocity field and the sym- 
metric positive definite matrix <x represents the normal- 
ized conformation tensor of polymer molecules and 1 is 



the unit tensor. The solvent viscosity is denoted by v and 
rj is the zero-shear contribution of polymers to the total 
solution viscosity v% = v(l + rj) and is proportional to the 
polymer concentration. In absence of flow, u = 0, poly- 
mers relax to the equilibrium configuration and <x = 1. 
The trace trcr is therefore a measure of polymer elonga- 
tion. 

The simplest geometrical setup that will prove use- 
ful to study the elastic turbulence regime for viscoelastic 
flows is the periodic Kolmogorov flow in two dimensions 
[8j. With the forcing / = (F cos(y / L) , 0), the system 
of equations |T][2|) has a laminar Kolmogorov fixed point 
given by 



u = (U cos(y/L),0) 
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with F = [i/U (l + ri)]/L 2 [§]. The laminar flow fixes a 
characteristic scale L, velocity Uo and time T = L/Uq. In 
terms of these variables, wc define the Reynolds number 
as Re — and the Weissenberg number as Wi = . 
The ratio of these numbers defines the elasticity of the 
flow El = Wi/Re. 

It is well known that the Kolmogorov flow displays 
instability with respect to large-scale perturbations, i.e. 
with wavelength much larger than L. In the Newto- 
nian case, the instability arises at Re c — y/2 [lfj]. At 
small Reynolds numbers, the presence of polymers can 
change the stability diagram of laminar flows [ll|, G3 
or induce elastic instabilities which are not present in 
Newtonian fluids [^, [l3|, [3] ■ The Kolmogorov flow is no 
exception, and recent analytical and numerical investi- 
gations have found the complete instability diagram in 
the Re-Wi plane @. For the purpose of the present 
work, we just have to recall that linear stability analysis 
shows that for sufficient large values of elasticity, the Kol- 
mogorov flow displays purely clastic instabilities, even at 
vanishing Reynolds number (see Fig. 1 of 0). We re- 
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mark that the fact that the original flow has rectilinear 
streamlines does not exclude the onset of the elastic insta- 
bility. Above the elastic instability the flow can develop 
a disordered secondary flow which persists in the limit of 
vanishing Reynolds number and eventually leads to the 
elastic turbulence regime [15j . 

The equations of motion (|H2p are integrated by means 
of a pseudo-spectral method implemented on a two- 
dimensional grid of size Lq — 2n with periodic bound- 
ary conditions at resolution 512 2 . Nmerical integrations 
of viscoelastic models are limited by Hadamard instabil- 
ities associated with the loss of positiveness of the con- 
formation tensor [l|| • These instabilities are particularly 
important at high elasticity and limit the possibility to 
investigate the elastic turbulent regime by direct imple- 
mentation of equations (|f 1 2j) . To overcome this problem, 
we have implemented an algorithm based on a Cholesky 
decomposition of the conformation matrix that ensures 
symmetry and positive definiteness (l7j and allows to 
reach high elasticities. 
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average power injection for the turbulent flow is simply 
P FL 2 
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FIG. 1: Mean power injection normalized with its laminar 
value r = P/Piam as a function of Wi = tU/L for a set of 
simulations with v = 0.769, r\ = 0.3, L = 1/4 and r = 4. 
The elasticity is El = 64 and the maximum Reynolds number 
Re — Wi/El is Re — 0.7. The line represents the power- 
law behavior r ~ Wi ' 25 . Inset: amplitude of the Reynolds 
stress II r (x), polymer stress H p (+) and viscous stress IR 
(*) nondimensionalized by the total stress amplitude FL. 

One of the main features of a turbulent-like regime is 
the growth of the flow resistance to external forcing. This 
can be quantified as the power needed to maintain a given 
mean velocity in the turbulent flow. The power injection 
in ([TJ) is P = (f • u) which, for the laminar flow © be- 
comes Piam — U^v{l + rj)/(2L 2 ). A remarkable feature of 
the Kolmogorov flow is that even in the turbulent regime 
the mean velocity and conformation tensor are accurately 
described by sinusoidal profiles [HI]: (u x ) = U cos(y/L), 
{<?xy) = — ^sin(y/L) with different amplitudes with re- 
spect to the laminar fixed point. Therefore the reduced 
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Figure [1] shows the behavior of the power injection as a 
function of the Weissenberg number Wi = tU/L. We see 
that at Wi ~ 15 there is a transition from the laminar 
regime to a turbulent-like regime in which r > 1. The 
growth for the higher values of Wi is compatible with a 
power law scaling r ~ Wi 25 which is qualitatively sim- 
ilar with experimental observations [l9l | where an expo- 
nent 0.49 is found. The different exponent observed here 
can be ascribed to the two-dimensional nature of our flow 
or to its geometrical property (rectilinear streamlines and 
absence of material boundaries). 

Because the Reynolds number in Fig[T]is always small, 
and therefore the inertial term in (fTJ) is negligible, it is 
natural to ask where is the origin of turbulent fluctua- 
tions. The momentum budget, in stationary conditions, 
reads 

dyU r = dy(n V + U p ) + f x (5) 

where IR = (u x u y ) is the usual Reynolds stress, IR = 
vd y {u x ) the viscous stress and Tl p = 2vrjT~ 1 (a xy ) is 
the stress induced by polymers. The numerical obser- 
vation that also the Reynolds stress is well described by 
a monochromatic profile, (u x u y ) = U2 sin(y/L), allows 
us to write the momentum budget for the amplitudes 
as FL = U 2 + vUjL + (2i/r?/r)S. The inset of FigQ] 
shows the different contributions (normalized with the 
total stress) as a function of Wi. In the laminar regime 
(Wi -> 0) U 2 = and from © one has Ilp/IR = 77. 
Above the transition to elastic turbulence, the polymer 
stress starts growing and reaches a value close to the vis- 
cous stress at the present maximum Weissenberg number. 
We remark that we observe no indication of saturation 
and therefore we may expect n p to become the domi- 
nant term at larger values of Wi. The contribution of 
the Reynolds stress always remains smaller than 10 ~ 2 , 
confirming the irrelevance of inertial terms. This is the 
hallmark of elastic turbulence where elastic stress has the 
role played by the Reynolds stress in usual turbulence. 

In order to get more insight in the elastic turbu- 
lence flow, in Fig. [2] we show two snapshots of the two- 
dimensional vorticity field at two different Wi. The first 
snapshot is taken at Wi = 22, slightly above the elas- 
tic instability threshold. The flow in this regime is still 
not turbulent and a secondary flow in the form of thin 
filaments is clearly observable. These small scale fila- 
ments, moving along the x direction, are elastic waves, 
reminiscent of the Alfven waves propagating in presence 
of a large scale magnetic field in plasma. Indeed, the 
possibility of observe elastic waves in polymer solution 
was theoretically predicted within a simplified uniaxial 
elastic model @ which has strong formal analogies with 
MHD equations, but they were never observed before. 



3 




FIG. 2: (color online) Snapshot of vorticity field at Wi = 
22 (left) and Wi = 45 (right). The flow is forced with a 
Kolmogorov forcing f x = Fcos(y/L) with L — 1/4. 



At higher values of elasticity the vorticity pattern be- 
comes progressively more irregular with chaotic motion 
of filaments. At Wi = 45 we observe a highly irregular 
pattern in which the underlying basic flow is hardly dis- 
tinguishable. This is the regime of elastic turbulence in 
which the flow develops active modes at all the scales. 
Fig. [3] shows the power spectrum of velocity fluctuations 
averaged over several configurations like the one shown in 
Fig. [21 A power-law behavior E(k) ~ k~ a is clearly ob- 
servable with a spectral exponent a larger than 3. Again, 
this is in quantitative agreement with what observed in 
laboratory experiments [ij and with the theoretical pre- 
dictions based on the uniaxial model [6j|. 
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FIG. 3: Velocity fluctuation spectra at Wi = 26 (squares) 
and Wi = 45 (circles). The line represents the power-law 
behavior k~ S 8 . 

One of the most promising applications of elastic tur- 
bulence is efficient mixing at very low Reynolds number. 
This is an issue of paramount importance in many in- 
dustrial problems, namely in microfluidic applications. 
Indeed, laboratory experiments in curvilinear channels 
have demonstrated that very viscous polymer solutions 
in the elastic turbulence regime are very efficient for small 



scale mixing Q. Mixing efficiency of polymer solutions 
has been studied in various setups, including microchan- 
nels 0] and two-dimensional magnetically driven flows 
Because in the elastic turbulent regime the flow is 
smooth (i.e. the energy spectrum is steeper that k~ 3 ) a 
suitable characterization of mixing is given in terms of 
Lagrangian Lyapunov exponent A [20J . This is defined as 
the mean rate of separation of two infinitcsimally close 
particles transported by the flow and, in the present case, 
is related to the polymer stretching rate [211 ]. 
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FIG. 4: Lagrangian Lyapunov exponent as a function of the 
Weissenberg number. Inset: Cramer function for three dif- 
ferent Weissenberg numbers: Wi = 12 (solid line), Wi — 21 
(dashed line) and Wi = 53 (dotted line). 

Figure 2] shows the behavior of the Lyapunov expo- 
nent as a function of Wi at fixed Re = 1. We observe 
that, above Wi — 10, A grows and saturates at values 
~ 2/r for Wi larger than 60. We remark that this be- 
havior is opposite to the one observed in the case of high 
Reynolds numbers viscoelastic flows where the injection 
of polymers reduces the degree of chaoticity by lowering 
A below 1/r [21 1. 

In the inset of Fig. 2] we plot the Cramer function G(^) 
which is defined from the probability density functions 
of finite-time Lyapunov exponents ft (7) ~ exp(— tG("f)) 
[ioj ]. As it is evident, increasing Wi not only the degree 
of mean chaoticity increases, but also fluctuations be- 
comes larger, in particular the distribution of 7 becomes 
asymmetric with a larger relative probability of positive 
fluctuations. It is remarkable that the same qualitative 
behavior is observed in the case of high- Reynolds New- 
tonian turbulence, where the distribution of Lyapunov 
fluctuations becomes more asymmetric with increasing 
Re Q. This suggests that in elastic turbulence elastic- 
ity (i.e. Wi) plays a similar role as non-linearity (i.e. Re) 
in ordinary hydrodynamic turbulence. 

Finally, we have investigated the dependence of poly- 
mer statistics on the Weissenberg number. In Figure [5] 
we show the average squared polymer elongation (trtr) 
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integrated over the flow volume and the amplitude of 
cross stress E. At small Wi these follow the laminar 
behavior, i.e. 2 + Wi 2 /4 and Wi/2, respectively. At 
the onset of elastic turbulence the cross polymer stress E 
grows faster than linearly in Wi, as already shown in Fig- 
ure [TJ and eventually appears to approach a power-law 
behavior with a slope close to 1.5. On the contrary, the 
squared polymer elongation in clastic turbulence grows 
more slowly than its laminar value. This is probably 
due to the loss of coherence in stretching experienced in 
the turbulent flow. At large Wi an asymptotic behavior 
~ Wi 15 appears to set in for tnx as well and the ratio 
(trcr) /E becomes constant. 
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FIG. 5: Average polymer elongation {trcr) (open circles) and 
polymer stress amplitude E (filled circles) as a function of Wi. 
Lines represent the laminar behaviors 2 + Wi 2 /4 (continuous) 
and Wi/2 (dotted). 

Summarizing, we have shown that elastic turbulence 
can be successfully reproduced numerically with the aid 
of a widely known viscoelastic model of polymer solutions 
(the Oldroyd-B model) and a simple geometrical setup 
(the two-dimensional Kolmogorov flow). Most observed 
features have a strong qualitative resemblance with ex- 
perimental results. Quantitative differences exist, how- 
ever, and may be traced back to the two-dimensional or 
to the boundaryless nature of our toy flow, or both. In 
this context it would prove extremely useful to perform 



numerical simulations in more realistic geometries and 
dimensionality to ascertain the origin of such differences. 
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